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We attempt to reconcile seemingly conflicting experimental results on the Higgs boson mass, the 
anomalous magnetic moment of the muon, null results in search for supersymmetry at the LHC 
within the 8 TeV data and results from _B-physics, all within the context of supersymmetric grand 
unifled theories. Specifically, we consider a supergravity grand unification model with non-universal 
gaugino masses where we take the SU(3)c gaugino field to be much heavier than the other gaugino 
and sfermion fields at the unification scale. This construction naturally leads to a large mass 
splitting between the slepton and squark masses, due to the mass splitting between the electroweak 
gauginos and the gluino. The heavy Higgs bosons and Higgsinos also follow the gluino toward large 
masses. We carry out a Bayesian Monte Carlo analysis of the parametric space and find that it can 
simultaneously explain the large Higgs mass, and the anomalous magnetic moment of the muon, 
while producing a negligible correction to the Standard Model prediction for Br (5° — ^ /x"'"/i~). We 
also find that the model leads to an excess in the Higgs diphoton decay rate. A brief discussion of 
the possibility of detection of the light particles is given. Also discussed are the implications of the 
model for dark matter. 



I. INTRODUCTION 

The CMS and ATLAS collaborations have discovered 
and measured [THn] the mass of a new boson which is 
most likely the Higgs boson [BHS] responsible for break- 
ing electroweak symmetry. In supersymmetry, one would 
identify this as the light CP-even Higgs boson [TDHIS], 
/i°. Both experiments agree that the mass is between 
125 and 126 GeV. It is quite remarkable that the ob- 
served Higgs boson mass lies close to the upper limit 
predicted in grand unified supergravity models [T7H^ 
which is roughly 130 GeV [M US [1IH21]- (For a recent 
review of Higgs and supersymmetry see [25].) Because 
the mass of the /i° boson in supersymmetry [^BH28, is 
less than that of the Z boson at the tree level, a large 
loop correction is necessary to match the measured value. 
The dominant one-loop Higgs self energy correction arises 
from its coupling to the top supermultiplet so that 
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where v — 246 GeV, Afs is the average stop mass, 
Xt = At—^J, cot P, fiis the Higgs mixing parameter and At 
is the trilinear coupling (both at the electroweak scale), 
and tan/3 = {H2) / {Hi) , where H2 gives mass to the up 
quarks while Hi gives mass to the down quarks and lep- 
tons. Since Ato^o has a logarithmic dependence of Afs, 
a sizable Am^o correction implies that the scale Ms is 
high, lying in the several TeV region. 

A high SUSY scale is also suggested by the ATLAS 
and CMS collaborations. So far, the LHC has delivered 
5.3 fb"^ and 23 fb"^ of integrated luminosity gS] at 7 TeV 



and 8 TeV respectively to both CMS and ATLAS. Analy- 
sis of large portions of this data in search of supersymme- 
try has only yielded null results, though it is important 
to note that the parametric exclusion limits provided are 
typically only on minimal or simplified models. When- 
ever one works with non-minimal models of supersym- 
metry, it is necessary to evaluate the signal efficiencies 
specific to one's model and determine the credible region. 
The null searches can be evaded obviously by just rais- 
ing the masses of the superpartners, and thereby raising 
the scale of SUSY, but it can also be done by producing 
mass hierarchies and mass splittings that are atypical in 
minimal models. 

The search for the rare decay _B° -^ /^^M^ ^-Iso has 
important implications for supersymmetry. The LHCb 
collaboration has recently observed [30] this rare de- 
cay, determining the branching ratio Br (B'^ — > ji^ ^i^^ = 
(3.2l^i['2) X 10^^, which is in excellent agreement with the 
Standard Model, and thus requires the supersymmetric 
contribution [5TH55J to this decay to be very small. This 
contribution is mediated by the neutral Higgs bosons and 
will involve a fiavor-changing scalar quark loop. (It is also 
sensitive to CP violation [5H I35j.) In the large tan/3 
limit, the branching ratio is approximately |36l I37j 
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where tb^ is the mean lifetime, fs, is the decay constant, 
and V^'^f is the effective CKM matrix element. The loop 
factors eg ^nd ey are given in terms of soft breaking pa- 
rameters of the 3rd generation m 



Q' 



rrij^, which are 



the masses of the left-handed squark, up-type squark, 
and down-type squark, as well as the gluino mass rUg, 
the strong coupling constant as, and the CP-odd Higgs 
mass niA'- 
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We note that the branching ratio given by Eq. ^ is 
suppressed by the factor {mt/niA)^ and so a large weak 
scale of SUSY which implies a large ttia , naturally leads 
to a small contribution to Br (B° — >• /i+/i^). Addition- 
ally, we see in Eq. ^ the factor (tan^/50)^, which im- 
plies that the SUSY contribution to B^ — >• n~^fi~ is fur- 
ther suppressed if tan/3 < 50. Together these effects also 
reduce the SUSY contribution [36l [38| to Sr (5 ^ Xsj) 
to negligible value. 

While the observation of a high Higgs boson mass, null 
results on the discovery of sparticles and the observation 
of no significant deviation in the B'^ — >■ /x''"/i~ branching 
ratio from the Standard Model result all appear to in- 
dicate a high scale for SUSY, the opposite is indicated 
by the Brookhaven experiment E821 [3SJ which measures 
'^/i — 5(5^ ~ 2) to deviate from the Standard Model pre- 
diction [3ni m] at the 3 a level. If this deviation is taken 
to arise from supersymmetry, then 
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The SUSY contribution [^2H15] arises from x^-^V and 
X°-/I loops. A rough estimate of the supersymmetric 
correction is 

Ja^ ~ sgn(M2M) (130 x 10'") P""*^""^ ) tan/3 , (7) 
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where MgusY is the SUSY scale. In order to obtain a 
SUSY correction of size indicated by Eq. ^ the masses 
of sparticles in the loops, i.e., the masses of x"^ , x°, Ji, 
and Vfj, must be only about a few hundred GeV. 

Another result which may be a signal of SUSY con- 
cerns the excess seen in the diphoton decay rate of the 
Higgs, which is above the Standard Model prediction. 
This excess is parametrized by the signal strength 
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and is reported as R^ 



1.6 ± 0.4 at CMS [3] and 



R^j = 1.8±0.5 at ATLAS g]. The excess is not statisti- 
cally conclusive and can easily be attributed to a simple 
fluctuation or to QCD uncertainties [50j. Still it is worth- 
while to consider how SUSY can contribute to this loop- 
induced decay (considering /i° in place of H). The excess 
in the diphoton rate has been discussed in a variety of 



models by various authors (see, e.g., [T3J[3T1|52] and the 
references therein). Within the MSSM, the largest con- 
tributions would arise via a r triangle, provided that its 
mass is not too high. (We discuss the calculation of R^^ 
in more detail in Section V A ) So, if the diphoton result 
is real, we have another indication of low scale SUSY. 

Assuming that the 3^ — 2 and the diphoton rate hold 
up, one has apparently conflicting results for the weak 
scale of SUSY. On the one hand, the high Higgs bo- 
son mass, null results on the observation of sparticles 
at the LHC, and the lack of any significant deviation in 
the Br (i?^ — > ^i^ jij branching ratio from the Standard 
Model prediction point to a high SUSY scale, i.e., a SUSY 
scale lying in the several TeV range. On the other hand, 
the 3 (T deviation in a^ and a fledgling excess in the dipho- 
ton decay of the Higgs boson decay point to a low SUSY 
scale lying in the sub-TeV range. These results taken 
together, point to a split scale SUSY with one scale gov- 
erning the colored sparticle masses and the heavy Higgs 
boson masses, and the other SUSY scale governing the 
uncolored sparticle masses. To generate this split scale 
SUSY, we construct in this work a supergravity grand 
unified model |17H19j by introducing non-universalities 
in the gaugino sector with the feature that the gaugino 
mass in the 5'C/(3)c sector is much larger than the other 
soft masses. In this model, radiative electroweak sym- 
metry breaking [53H55J (for a review see |56j ) is driven 
by the gluino mass. In this work, we label this model 
as gSUGRA. We wih show that jSUGRA satisfies aU 
of the experimental results simultaneously by exploiting 
a feature of the renormalization group equations which 
leads to a splitting between the squarks, gluino, Higgs 
bosons, and Higgsinos which become very heavy, and the 
sleptons, bino and winos which are allowed to remain 
light at the electroweak scale. (The sfermion masses still 
unify at a high scale.) We will use a Bayesian Monte 
Carlo analysis of gSUGRA to show that it satisfies all 
experimental results and determine the credible regions 
in the parameters and sparticle masses. 

The outline of the rest of the paper is as follows: In Sec- 
tionJHJ we discuss the general framework of non-universal 
SUGRA models with specific focus on gSUGRA where 
the gaugino mass in the SU{?i)c color sector is much 



larger than other mass scales in the model. In Section HI 



we discuss the statistical framework used in our Bayesian 
Monte Carlo analysis of a simplified parametric space for 
gSUGRA. In Section IV we explore the impact of LHC 



searches for sparticles on gSUGRA using event-level data 
and signal simulations. The results of our analyses as well 
as the details of Higgs diphoton rate are presented in Sec- 
tion Iv] Concluding remarks are given in Section VI 



II. THE sSUGRA MODEL 



Supergravity grand unification [T7HT9] is a broad 
framework which depends on three arbitrary functions: 
the superpotential, the Kahler potential, and the gauge 



kinetic energy function. Simplifying assumptions on the 
Kahler potential and the gauge kinetic energy function 
lead to universal boundary conditions for the soft pa- 
rameters which is the basis of the model referred to as 
mSUGRA/CMSSM. The parameter space of mSUGRA 
is given by toq, m.1/2, ^o, tan/3, and sgn(/i), where 
Too is the universal scalar mass, mi/2 is the universal 
gaugino mass, Aq is the universal trilinear coupling, and 
tan/3 = {H2)/{Hi). Here H2 gives mass to the up quarks 
and Hi gives mass to the down quarks and the leptons, 
and /i is the Higgs mixing parameter which enters in the 
superpotential as /ii?ii?2- 

However, the supergravity grand unification frame- 
work does allow for non-universalities of the soft pa- 
rameters, i.e., non-universalities for the scalar masses, 
for the trilinear couplings and for the gaugino masseaM 
In gSUGRA, we consider supergravity grand unification 
with universal boundary conditions in all sectors except 
in the gaugino sector. In this sector, we specify that the 
SU(3)c gaugino mass, M3, be much larger than the uni- 
versal scalar mass and also much larger than the gaugino 
masses M2,Mi in the SU(2)i, U(l)y sectors, i.e., 

M3>mo,Mi,Af2 (9) 

The constraints of Eq. ^ ensure that the radiative 
breaking of electroweak symmetry will be driven by the 
gluino (hence, jSUGRA). Now, the gluino mass enters 
in the renormalization group equations for the squark 
masses and thus the squark masses will be driven to val- 
ues proportional to the gluino mass as we move down 
from the GUT scale toward the electroweak scale. Con- 
sequently, a gluino mass in the ten TeV region will also 
generate a squark mass in the several TeV region. On the 
other hand, the RGEs for the sleptons do not depend on 
the gluino mass at the one-loop level and if mo,Mi,M2 
are O(lOOGeV), the masses of the sleptons as well as 
the electroweak gauginos at the electroweak scale will 
likely remain this size. Thus the RG evolution creates 
a natural splitting of masses between the squarks and 
the sleptons at the electroweak scale even though they 
have a common mass at the grand unification scale. The 
renormalization of these soft masses for a sample point 
in gSUGRA is shown in Fig. IT] The huge mass split- 
ting between the squark and slepton masses at low scales 
even though they are unified at high scales is reminiscent 
of the gauge coupling unification where the three gauge 
couplings ai which are split at the electroweak scale but 
come together at the grand unification scale. We note 
that the split spectrum of gSUGRA is very different in 
nature from that of what is commonly called "split super- 
symmetry" [71], which consists of light Higgsinos Hu^d, 
B, W, g and one Higgs doublet but does not allow for 
light sfermions. 



In GUT models, non-universal gaugino masses can 
arise from superfields that transform as a non-singlet 
IRs of the GUT group and get VEVs in the spontaneous 
breaking and give masses to the gauginos. The general 
form of the gaugino mass term in the Lagrangian is 
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where {F)ab is a non-zero VEV of mass dimension 2, 
and Mpi is the Planck mass. The A's belong to the 
adjoint of the GUT group: 24 for SU(5) and 45 for 
SO(IO). Now only the symmetric product of the ad- 
joints enters in the analysis. Thus for SU(5) one has 

75 © 200, while for SO(IO) 
© 54 ® 210 ® 770. With the 



(24 (E) 24)sym = 1 ® 24 i 



^ The literature on non-universalities in SUGRA models is enor- 
mous. For a sample of early and later works see |57H69| and for 
a review see |70| . 



one has (45 (g) 45)syni = 
use of singlet and non-singlet breaking, one can produce 
a hierarchy in the gaugino masses so that Eq. (l9| holds. 
We note that non-universalities of gaugino masses arise 
also in string based models, see, e.g., [7^ . 

In our study of gSUGRA, we introduce gaugino sector 
non-universalities by having mi/2 — > rhi/2 = Mi = A/2 
and M3 = 10 mi /2 as an illustrative example, so that at 
the unification scale. Mi : M2 : M3 = 1 : 1 : 10. We now 
show how this choice can be constructed by combining 
singlet and non-singlet breaking in SU(5) and in SO(IO). 
In SU(5) we consider the linear combination 1 -I- a 24 -|- 
6 75. Now the singlet breaking gives the ratio Mi : M2 
A/3 = 1 : 1 : 1, the 24-plet gives the ratio [H] (-1/2 
—3/2 : 1) while the 75-plet gives the ratio [66] (—5 
3 : 1). Choosing a = —8/11 and h ~ —1/11 leads to 
the desired ratio A/i : A/2 : A/3 = 1 : 1 : 10. This 
scheme also applies to SO(IO) since SU(5) c SO(IO). 
However, for SO (10) we can also consider gaugino mass 
terms in representations of SU(4) x SU(2)i x SU(2)fl C 
SO(IO) and label the breaking terms by SU(4) x ^\]{2)r 
representations as subscripts. In this case we consider the 
breaking 1 + a210(i^i) -I- 5210(i5_i) where the 210(i_i) 
gives the gaugino mass ratio [BB] of (—3/5 : 1 : 0) and 
210(i5_i) gives the gaugino mass ratio pBl of (—4/5 : : 
1). Thus we can choose a = —3/4 and b = 3/2 to get the 
desired 1 : 1 : 10 ratio. We limit ourselves to this ratio 
for the rest of the analysis in this paper. However, many 
features of this analysis will persist with different ratios 
of A/i : A/2 : A/3 as long as A/3 > mo, A/i, A/2. 

In gSUGRA, radiative electroweak symmetry breaking 
is dominated by the large gluino mass which is responsi- 
ble for giving large masses to the squarks. We contrast 
this work with other recent works which have attempted 
to explain g^ — 2 in the context of a high Higgs boson 
mass. This is attempted in |51| with the assumption of 
a light slepton and heavy squark spectrum. The analy- 
sis also tries to correlate (7^ — 2 with the diphoton rate. 
However, this model is not a high scale model and the 
analysis is limited to assumptions of the spectrum at the 
electroweak scale. In [73] the authors assumed a split 
family supersymmetry. The analysis of 174] uses non- 
universal gaugino masses in an SU(5) model but the de- 
tails of the model are significantly different from the work 
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FIG. 1. Two-loop renormalization group evolution of the 
soft parameters in jSUGRA. The input parameters used here 
are those of the best-fit point determined from our analysis 
in Section |III| The fields are labeled in the figure and also 
in color. The gaugino fields are presented in black and the 
Higgs fields are presented in green. The squarks and sleptons 
are in blue and red, where the left-handed fields are solid and 
the right-handed fields are dot-dashed. Additionally, mo is 
the soft mass for the scalars, mi/2 is the common mass of the 
U(l)y and SU(2)_t gaugino fields, and ji is the Higgs mixing 
parameter. 



presented here. The work [TFj also addresses the issue of 
getting Hght uncolored and heavy colored particles but 
the analysis is within a gauge mediated supersymmetry 
breaking. 

The attractive feature of jSUGRA is that the rela- 
tively large value of M3 automatically drives the squarks 
to be massive while the sleptons as well as the bino and 
the light wino are left alone. This is illustrated in Fig. [l] 
where we display the renormalization group flow for a 
sample point from our analysis. We wish to show that 
this simple feature automatically satisfies all of the em- 
pirical results that we have discussed here that hint at 
the supersymmetric spectrum. To this end, we perform 
a Bayesian Monte Carlo analysis of jSUGRA with the il- 
lustrative example of the (1:1: 10) gaugino mass ratio, 
which we discuss in the sections that follow. 



III. STATISTICAL FRAMEWORK 

We study here the parameter space of gSUGRA for the 
case where the ratio of the gaugino masses at the GUT 
scale is 1 : 1 : 10. In this case, jSUGRA is parametrized 
by Too, TOi/2, Aq, and tan/3 (having selected sgn {jjl) = 1). 
Here, TO1/2 = Mi = M 2 while Afg = 10toi/2. The di- 
mensionful parameters toq, ^i/2: and ^0 are all speci- 
fied at the GUT scale. The ratio of the two Higgs VEVs 
tan/3 — {H2) / {Hi) , is specified at Mz- We further in- 
clude four Standard Model nuisance parameters to create 
an 8D parameter space. Namely, we add the top quark 



pole mass, the running bottom quark mass, the strong 
coupling, and the EM coupling. We create from these 
the parameter space 0: 



= |too,TOi/2, Ao,tan/3,TO^ 
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For each parameter Oi € 0, we begin by selecting uniform 
distributions in the allowed ranges prior to considering 
the experimental data. The prior distributions that we 
have selected for our parameters are uniform on either a 
linear or a log scale: 



mo e [50, 5000] GeV (log) 

TOi/2 e [50, 2500] GeV (log) 

Ao e [-50, 50] TeV (hnear) 

tan/3e [3,60] (linear) . 



(12) 



The nuisance parameters in are uniform in a 2 cr range 
(linear scale) around the central values, which are speci- 
fied in Table H 

Next we collect the relevant observables into D, which 
is a set of pairs of central values and uncertainties of ex- 
perimental measurements. The observables include the 
precise measurements of the nuisance parameters, along 
with the results from flavor physics Br (B^ — > /i"*'/i~) and 
Br {B -^ Xsj), the muon anomalous magnetic moment 
(5ap , the measured mass of the (ostensibly) light CP-even 
Higgs boson, as well as limits on superpartner masses. 
We further include the fit to the thermal relic density of 
dark matter, w^ = O^^/i^, from CMB temperature fluc- 
tuations measured by WMAP (9 year dataset) (SU] and 
Planck (15.5 month dataset) ^. In jSUGRA, the light- 
est neutralino is indeed a candidate for cold dark matter, 
but we wish to allow for multicomponent models of dark 
matter, and so we only consider the upper limit of lu^. 
The central values and uncertainties of D are specified in 
Table m 

The goal now is to update our a priori guess for the 
probability distributions of the parameters in (given in 
Eq. ( 12 1) with the empirical information in D, giving the 
posterior probability distribution. This distribution can 
then be marginalized to determine the credible region of 
one or two parameters. The calculation of the posterior 
probability distribution is achieved using Bayesian infer- 
ence, but we first need to be able to compare a parametric 
point in our model to the empirical data in D. This re- 
quires a set of mappings ^^ : — >■ M corresponding to 
each di e D, which just give the theoretical calculation 
for the observable corresponding to each di. These map- 
pings are computed using numerical codes incorporated 
in our analysis software SusyKit [52] . 

Now we can move on to constructing the pos- 
terior probability distribution, which is given by 
Bayes' theorem 



P(0|D) 



P(D|0)P(0) 
P(D) 
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Observable 


Central value 


Exp. Error Th. Error 


Distribution 


Ref. 


SM Nuisance Parameters 


mr'" (GeV) 


173.5 


1.0 


Gaussian 


m 


mtimt)^ (GeV) 


4.18 


0.03 


Gaussian 


m 


as(Mz)^ 


0.1184 


7x10"* 


Gaussian 


m 


„ -1 r tf ^MS 


127.933 


0.014 


Gaussian 


m 


Measured 


5a^ X 10" 


287 


80 10 


Gaussian 


[39H1I] 


Br {B° ^ ^i+ fi-) X 10» 


3.2 


1.92 14% 


Gaussian 


ED] 


Br [B ^ Xs-y) X 10^ 


3.55 


0.26 0.21 


Gaussian 


m 


Br {B+ -^ T+v) X 10'' 


1.79 


0.48 0.38 


Gaussian 


m 


^x 


0.1126 


0.0036 10% 


Upper-Gaussian 


m 


/i° Mass (GeV) 


125.7 


0.2 2.0 


Gaussian 


mm 


95% CL Particle Mass Limits (GeV) 


/i" 


122.5 


- 


Lower - Step Func. 


m 


h° 


129 


- 


Upper - Step Func. 


m 


xS 


46 


5% 


Lower - Error Func. 


m 


xl 


62.4 


5% 


Lower - Error Func. 


m 


xi 


99.9 


5% 


Lower - Error Func. 


m 


x2 


116 


5% 


Lower - Error Func. 


m 


xf 


94 


5% 


Lower - Error Func. 


m 


BK 


107 


5% 


Lower - Error Func. 


m 


m 


94 


5% 


Lower - Error Func. 


m 


n 


81.9 


5% 


Lower - Error Func. 


m 


6i 


89 


5% 


Lower - Error Func. 


m 


ii 


95.7 


5% 


Lower - Error Func. 


m 


5 


500 


5% 


Lower - Error Func. 


m 


9 


1100 


5% 


Lower - Error Func. 


m 



TABLE L Summary of the observables used to construct the likelihood function. The distribution labeled "Upper-Gaussian" 
used for the oj^ observable means that there is only a decrease in likelihood for values larger than the central value. The 
95% CL limits are evaluated using the complementary error function, as the bound is smeared by the theoretical uncertainty. 
Limits specified with a step function distribution indicate a hard cut, where points on the wrong side of the limit are assigned 
zero likelihood. 



P(0) is the prior distribution given in Eq. (12). The de- 
nominator is the so-called Bayesian evidence Z = P{D), 
which can be used in model selection tests, but as we 
are only interested in parameter estimation, it serves as 
a normalization constant. The final factor is the likeli- 
hood function £ — P(D|0), which is constructed by the 
"pulls" method 



21n£ 



E 



(6(0) 



(14) 



where <Ti and t^ are the experimental and theoretical un- 
certainties, respectively. This is straightforward for the 
case that a measurement with precision is reported. In 
many cases only the 95% CL limits are given. In those 
cases, a smearing due to the implicit theoretical uncer- 
tainty in the computation is used and the likelihood is 
computed from the complementary error function. A 
hard cut on an observable can also be made by using 
a step function, i.e. assigning zero likelihood to points 
that are on the wrong side of a limit. The numerical val- 
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FIG. 2. A display of the mass spectrum for sparticles and the Higgs boson with split scales, i.e., a low scale for Xii Xi i "^i, ■?2, ^ 
and a high scale for H^ ,X2 i^iiQiS- Shown are the credible intervals in the superpartner masses from the Bayesian analysis 
of jSUGRA. The lighter superpartners are presented in the left panel, and the heavier are presented in the right panel. The 
posterior means are indicated in red. 



lies used to construct the likelihood function is given in 
Table n 

Our analysis was performed using our software pack- 
age SusyKit [S2], which uses the efhcient multi-modal 
ellipsoidal nested sampling algorithm implemented in 
the MultiNest [83-85^ Ubrary. Additionally, SusyKit 
interfaces with several standard numerical codes such 
as SOFTSUSY [SS], MicrOMEGAs WT, 'SSi, Feyn- 
HiGGS [521 [10], and SuperIso Relic [91, 92 . SusyKit 
is written entirely in C++11 and is largely inspired by 
the FORTRAN-90 code SuperBayes [MllM]. 

We specify the MultiNest sampling parameters 
niive = 5,000 and tol = 0.01. The analysis has required 
the evaluation of the likelihood function at 1.1 million 
points to sufficiently explore the parametric space. The 
result is a chain of 81,000 Monte Carlo sample points 
which is used to compute ID and 2D marginalized distri- 
butions in our principal and derived parameters, and to 
establish credible regions in these parameters. We found 
that the credible regions entered areas that would be ex- 
cluded by the LHC in minimal SUSY GUT models such 
as mSUGRA, so we found it necessary to evaluate the 
impact of LHC searches on gSUGRA. 



IV. LHC ANALYSIS 



In order to evaluate the impact of null results in the 
searches for supersymmetry at the LHC on gSUGRA, we 
construct an auxiliary likelihood function, £lhC) based 
on the Monte Carlo event generation and detector simu- 
lation for our sample points. 

We begin by generating 200,000 events for each sample 
point in our chain using PYTHIA [SS] [55] considering 



2 — > 2 SUSY production processes with y^s" = 8 TeV. 
We find that the total cross section for these processes is 
O (100 fb) and the dominant modes involve the produc- 
tion of X]*; X2; Xi 7 ^^ '^1^ ^2, and, £ and v. This is to 
be expected because in gSUGRA, the scalar quark fields 
all become heavy as they are renormalized to the elec- 
troweak scale, while the scalar leptons are allowed to re- 
main light to produce contributions to 6a^ and the Higgs 
diphoton decay rate. By investigating the dominant de- 
cays of these particles, we decide that supersymmetry 
searches in leptonic final states are the most relevant to 
gSUGRA. We have used the 3£ and same-sign 2i searches 
at CMS [97] using 9.2 fb"^ at ^^ = 8 TeV to construct 
our £lhc- These searches are performed using 108 and 
4 event bins respectively, which serve as counting exper- 
iments and are naturally Poisson distributed. Therefore 
-Clhc is computed by 



-Clhc ^ 11 -^i 



(15) 



iGbins 



Each £,i would be a simple Poisson likelihood, except 
that one of the parameters to the Poisson distribution, 
the expected background yield, bi , can have a large uncer- 
tainty, 6bi. Thus, it is necessary to convolve the Poisson 
distribution with a distribution for the background yield. 
Naively this would be a Gaussian distribution, however in 
the case that the relative error in the background yield 
is large, i.e., Sbi/bi > 20%, then a non-trivial portion 
of the convolution is due to contributions from negative 
bi, or even if the integration is limited to non-negative 
background, a large portion of the PDF may be omitted. 
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FIG. 3. An exhibition of the 1 ct and 2 a credible regions of the marginalized posterior probability distributions for the 
parameters of interest of gSUGRA. Left panel: the credible regions in mo and rhi/2- Right panel: the credible regions in the 
dimensionless parameter Ao/mo and tan/3. The location of the best-fit point is indicated by a circled 'X' and the posterior 
mean is given with a solid dot. 
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FIG. 4. A display of the marginalized posterior probability distributions for gSUGRA in the parameters of interest as well as 
some important derived quantities. The top row (left to right) gives the posterior PDF for mo, fhi/2, ^o, and tan/3, and the 
bottom row (left to right) displays the same for the top quark mass, the light CP even Higgs boson mass, the contribution to 
— 21n£ due to the anomalous magnetic moment of the muon (which we have denoted as •x^{5a^), and the thermal relic density 
of cold dark matter, ujy^. The location of the best-fit point is indicated by a circled 'X' and the posterior mean is given with a 
solid dot. 



Thus as a heuristic, we use the following definition for 2,i 



Pois(si -I- b\ Oi)F{bi,5bi;b) 



(16) 



where i is the event bin, Pois is the Poisson probability 
mass function, Si is the expected signal yield, Oi is the 
number of observed events, and as defined already bi is 
the expected background yield, and 5bi is the uncertainty 



in the background. The function F is defined according 
to our heuristic 



F{b,,6h-b) = 



\Af{h,5h]b) , 5b,/bi <20% 
\\nN{b,,5h]b) ,6bi/b^> 20% 



(17) 



where Af is the Gaussian distribution and In A/" is the log- 
normal distribution. As a further heuristic, it is necessary 
to account for cases when either 6^ = or 6bi = 0. These 
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which is the contribution to — 21n£ due to (5a^. The level 
curves are given in the plane of the top mass and hP mass. The 
level curves are constructed by interpolating equally-weighted 
sample points. 



cases are clearly oversights in the CMS preliminary anal- 
ysis summary; still they must be addressed. We choose 
a sentinel value A = 10~^ and use Shi = A if Sbi is zero 
and we set bi = Sbi if hi is zero. 

The expected signal yield Si is the product of the ef- 
ficiency Ci with the total SUSY cross section and the 
integrated luminosity. The efficiency e^ is the propor- 
tion of the total generated events that would be counted 
in the i bin, and is determined by running the events 
through a detector simulation, which we have carried out 
with PGS4 IHSj. Jet objects were reconstructed using 
the anti-fcx algorithm, with a distance parameter of 0.5. 
We implemented the cuts to place events into bins in a 
modified version of Parvicursor [HS]. The object se- 
lection criteria, event vetoes, and geometrical cuts are 
reproduced as in [97] . 

To combine the likelihood from these searches to the 
likelihood function described in Section III we first com- 



pute the likelihood for the Standard Model according to 
this analysis by turning off the signal, £sm — •ClhcIs^q- 
We then add the likelihood ratio statistic to the full like- 
lihood function. 



21n£ 



-21n£ 



2 In min 



J -Clhc 

i -CsM 



(18) 



which is approximately x^ distributed, and is a natu- 
ral addition to the other "pull" terms in our likelihood 
function. Having computed the updated likelihood due 
to these CMS searches, it is necessary to re- weight the 
samples by a factor exp (Aln£). We can now proceed to 
determine the marginalized posterior probability distri- 
butions within our parameters of interest. 



V. RESULTS 

In this section we present the results from our 
Bayesian analysis. Given our likelihood function, we 
determine the Bayesian evidence of jSUGRA to be 
InZ = —11.9 ± 0.042. We provide this for reference, 
as we do not perform a model selection test. The 
best-fit point in our analysis is determined to have 
2.73, and leaving out some of the nuisance pa- 



Amin 



polc\ 



rameters, is specified by (77io,mi/2,^Oj tan/3,r7i( ) 
(341,429,298,9.73,174) where the massive parameters 
are specified in GeV. This point illustrates the gen- 
eral result of gSUGRA that high /i° mass and (Ja^ 
can be simultaneously satisfied. Additionally, the large 
scalar quark and gluino masses allow for consistency with 
Br (S° -^ M+/i") and Br {B -^ X,j). The credible re- 
gions in the masses of the heavier particles in jSUGRA 
are presented in the right panel of Fig. [2| and the light 
particles of ^SUGRA that create the tJa^ contribution as 
well as the contribution to the diphoton Higgs decay are 
given in the left panel. 

The 1 a and 2 a credible regions in our parameters of 
interest are given in Fig. [3] where we have chosen to use 
the dimensionless parameter Ao/mQ. The ID posterior 
distributions in these parameters are given in the top 
panels of Fig. |4j though here we did give the distribution 
for the dimensionful parameter Aq. 

While gSUGRA largely achieves the correct h^ mass 
and dttf^ contribution as shown in the middle two lower 
panels of Fig. |4J the posterior distribution in the top 
mass is shifted up from the central value by 0.5 GeV to 
174 GeV, which is evident in the lower left panel of Fig. [4] 



The tension between the top mass, the h^ mass and 5a^ 
is clearly displayed in Fig. [6] where we have interpolated 
sample points from a slice in our likelihood function and 
presented level curves in "x^(^a/i)" which is the contribu- 
tion to —2 In £ due to (5a^. It is evident that the higher h'^ 
mass and (5a^ is best matched in gSUGRA for a slightly 
heavier top quark. 

We point out that this tension is not overly significant 
in gSUGRA for two reasons. First, there is a large theo- 
retical uncertainty in the calculation of the /i° mass at the 
2-loop level, which when considered does lift most of the 
tension. Next, we specified in jSUGRA M3 = lOmi/2, 
where 10 is an arbitrary choice. Allowing the coefficient 
to be a new degree of freedom or simply selecting several 
different choices will likely resolve this tension as well. 

In our Bayesian analysis, we have sampled the param- 
eter space using the older WMAP7 value for uj^ in £ but 
we can see from the fourth panel from the left in the bot- 
tom row of Fig. |4] that the slightly larger value indicated 
by WMAP9 and Planck would simply enlarge our credi- 
ble region. Additionally, we see in Fig. [5] that gSUGRA 
is not currently constrained by the best available limit on 
the direct detection of xl dark matter, and is slightly be- 
yond the projected sensitivity of XENONIT and Super- 
CDMSIT, creating a sort of nightmare scenario for dark 
matter experiments, as our dark matter signal would be 
competing with the cosmic neutrino background. The 
LSP in our model is consistently a bino, and the x§ is 
a wino. There is virtually no mixing with the Higgsino 
sector as the Higgsino mass parameter fi becomes very 
large due to the large M3. The sensitivity to dark mat- 
ter experiments can be increased by adjusting the ratio 
of Ml to M2 to allow for greater bino-wino mixing within 
the LSP state. 

One of the exceptional aspects of gSUGRA is the 
presence of many light superpartners that have thus 
far evaded detection at the LHC. We concede that the 
searches that we considered here are not by any means 
comprehensive, but they are designed to constrain the 
production modes most prevalent in jSUGRA. The lim- 
its are evaded largely due to the stringent selection crite- 
ria and the difficulty in identifying r leptons. Addition- 
ally, the mass hierarchy of jSUGRA limits the possibility 
of cascading decays. 

We note that the parametric space of gSUGRA, natu- 
rally fits into the Hyperbolic Branch |100H102j of radia- 
tive breaking of the electroweak symmetry. This is due 
to the fact that the stop masses are driven to be large 
by the giuino, giving a large Q — JTfi^ mj , and it was 
shown in |1031 I104J that Q > 1 TeV corresponds to a hy- 
perbolic geometry of soft parameters that give radiative 
EWSB (a large SUSY scale in the tens of TeV also arises 
in a certain class of string motivated models |1051I106] ). 
Still, gSUGRA as it stands produces a large value of ^ 
with respect to the Z mass. Specifically, a large value 
of ^ is necessary to balance the large value of AI^ which 
enters in the corrections to the H2 field mass. 



A. Higgs Diphoton Decay 

In the Standard Model, the loop-induced decay of the 
Higgs into two photons is mediated mainly by the W, 
top, and to a lesser extent, the bottom quark. The partial 
width reads [TUT] 



r (ij -^ 77) = 



"'EM'"' if 

256w27r3 



J2 N,jQ}A,/2iy) + Ai{Tw) 

f=t,b 

(19) 
where Ti = Amf/m'jj, and the spin form factors are 

A/2(r)-2r(l-(r-l)/(r)) (20) 

^i(r) = -(2-f3T-3T(r-2)/(T)) (21) 

and the universal scaling function /(r) is 



arcsm r 



-1/2A 



f{r) = 



4 V i-^n^ 



r > 1 



r < 1 



(22) 

Supersymmetry corrects this partial width [28] by factors 
involving the Higgs mixing angle a and /3 arising from 
the two Higgs doublets. Additionally, new amplitudes 
are available mediated by the charged Higgs, charginos, 
and sfermions. The couplings to the charginos arise from 
Higgsino-gaugino mixing, but in gSUGRA the Higgsinos 
are very heavy thus the lighter chargino is always purely 
charged wino while the heavier one is purely charged Hig- 
gsino. This means that overall the chargino contribution 
is small either because the coupling is suppressed or be- 
cause the mass is too large. The charged Higgs exchange 
is also suppressed due to its large mass. Thus the largest 
contributions can come only from the sfermion sector, 
which in grSUGRA is dominated by the staus. 

In the decoupling limit where Ma ^ Mz which corre- 
sponds to a = P — 7r/2, the Higgs coupling to the staus 
is given by [3T1 [52] 



gii^TiTi — I^Ci =F Qt sin 6'w cos 29^ 
mT~{Ai — /i tan/3) 



T 



2Af| 



sin 2..-^ 



(23) 



with ci — cos^0y, and C2 — sin 9^. The '— ' case corre- 
sponds to i = 1, and the '-f' case corresponds to i = 2. 
The partial width in gSUGRA including the amplitude 
due to staus then reads 



r {h' ^ 77) = 



256w27r3 



J2 NcjQ}A,/,iTj) 



f=t,b 



Mi 



-Ai{tw)+ ^ 9h0riT2 f-^o(Ti) 



i=l,2 



(24) 



10 




I 1 (T Credible Interval 
I 2 a Credible Interval 



1 



1.4 



1.5 



1.6 



FIG. 7. A display of the marginalized posterior probability 
density of i?^.^ from our analysis. The 1 a and 2 a credible 
intervals are indicated in darker blues. We define R-y^ as the 
ratio of the diphoton partial width of the light CP-even Higgs 
boson to the corresponding width for a Standard Model Higgs 
of the same mass (see Eq. ([8|). 



and the spin zero form factor is 

Ao(r) = -r(l-r/(r)) 



(25) 



We identify the ratio of this partial width to the Standard 
Model width given in Eq. (19 1 as R^^. (We have taken 
the ratio of the theoretical and observed /i° production to 
be unity.) We compute this ratio for each of our Monte 
Carlo samples and construct the ID posterior PDF in 
this derived parameter which we present in Fig. [7] We 
find that gSUGRA generically produces a ~ 20% boost to 
this decay mode over the Standard Model case. The 2 a 
credible interval is [1.03, 1.38], which is quite consistent 
with the preliminary results arriving from the LHC. 



VI. CONCLUSION 

The recent observation of the Higgs boson mass around 
125 GeV points to large loop corrections which can be 
achieved with a large weak scale of SUSY. A large SUSY 
scale also explains the suppression of SUSY contributions 
to the decay B^ — >■ fi^n~, to be consistent with the re- 
cently measured branching ratio for this process. On the 
other hand, the experimental observation of a 3a effect 



in Sa^ and a possible excess in the diphoton rate R^j in 
the Higgs boson decay over the standard model predic- 
tion cannot be explained with a high SUSY scale. Thus 
the two sets of data point to a two scale SUSY spec- 
trum, one a high scale consisting of colored particles, 
i.e., the squarks and the gluinos, and the Higgs bosons 
(aside from the lightest Higgs) and the other a low scale 



for masses of uncolored particles including sleptons and 
the electroweak gauginos. 

In this work we discuss the high scale supergravity 
grand unified model, gSUGRA, which includes the fea- 
ture of a two scale sparticle spectrum where the spar- 
ticle spectrum is widely split at the electroweak scale. 
This is accomplished within supergravity grand unifica- 
tion with non-universal gaugino masses such that M3 3> 
Mi,M2,mQ. As an illustration we consider the specific 
case where Mi : M2 : -M3 = 1 : 1 : 10 at the unifica- 
tion scale, Ml = M^ = fhxi^ and M-^ >> mg. This case 
is designed to be mainly illustrative and can be easily 
embedded within SU(5) and SO(IO). Using a Bayesian 
Monte Carlo analysis. It is found that this construction 
simultaneously explains the high /i*^ mass, null results 
for squarks and gluino searches at the LHC, a negligible 
correction to the branching ratio for B^ — >■ /i^/i^ , a 3ct 
deviation of 5^ — 2 from the Standard Model prediction 
as well as the nascent excess in the diphoton signal of the 
Higgs. 

The observable sparticle spectrum at the LHC in this 
model consists of light sleptons and light electroweak 
gauginos. However, sleptons and electroweak gauginos 
are typically difficult to observer at the LHC and thus 
far have evaded detection in multi-lepton searches in ex- 
periments at the ATLAS and the CMS detectors with 
the 7 TeV and 8 TeV data. The most promising 2 -> 2 
processes that can generate sparticles at the LHC in this 
model are pp -^ Xi Xi^j X2X1 ■ The identifying signatures 
of such processes will indeed be multi-leptons and miss- 
ing energy. It is hoped that at increased energies and 
with larger luminosities such signals will lie in the observ- 
able region. However, a detailed analysis of the signals 
in needed requiring a knowledge of the backgrounds for 
these processes. 

Another aspect of the simplified jSUGRA model re- 
lates to the spin-independent x° — p cross section. This 
cross section is found to be rather small for the case when 
the gaugino masses are chosen in the ratio (1:1: 10). 
The reason for this smallness is easily understood. The 
constraint Mi = M2 at the GUT scale, leads to an 
LSP which is essentially purely bino with very little 
Higgsino or wino content. The purely bino nature of 
the LSP leads to a suppressed Xi — p cross section (see 
e.g., |108| ) which lies beyond the reach of the current 
and projected sensitivities for direct-detection experi- 
ments. However, the above result is very specific to the 
Ml : M2 : M3 = 1 : 1 : 10 assumption and a modification 
of the above should allow xj — p cross section within the 
observable range in the projected sensitivities for direct- 
detection experiments. We note that while our analysis 
was performed using the older WMAP7 measurement of 
the cold dark matter relic density, the newer measure- 
ments from WMAP9 and Planck (with 15.5 months of 
data) only slightly increase the measurement. As we only 
apply the upper limit from these measurements to allow 
for the possibility of multi-component theories of dark 
matter, the newer results would only expand the credible 



11 



regions of our parameter space and either increase or not 
affect at all the likelihood of our best-fit point. 



Finally, we note that the large squark masses in 
gSUGRA would also help stabilize the proton against de- 
cay from baryon and lepton number violating dimension 
five operators [1041 11091 1110) (for a review see |111) ). 
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